Pseudogap due to Antiferromagnetic Fluctuations and the Phase Diagram of 
High- Temperature Oxide Superconductors 
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A reduction of the density of states near the Fermi energy in the normal state (pseudogap) of high- 
temperature oxide superconductors is examined on the basis of the two-dimensional tight-binding 
model with effective interactions due to antiferromagnetic fluctuations. By using antiferromagnetic 
correlation lengths which are phenomenologically assumed, the doping dependence of the pseudogap 
is obtained. The superconducting transition temperature decreases and eventually vanishes due to 
the pseudogap as the hole concentration is reduced. 



PACS numbers: 74, 74.25-q, 74.25.Dw, 74.62.-c 

In high-temperature oxide superconductors (HTSCs), 
behaviors which can be attributed to a reduction of the 
density of states (DOS) near the Fermi energy (pseu- 
dogap) have been observed. For example, experimen- 
tal results in photoemission spectroscopies ||^~^, tunnel- 
ing spectroscopies 
NMR experiments 



|7|, a specific heat measurement 

and neutron scattering [h4|. 
However, the origin of the pseudogap remains controver- 
sial: preformed Cooper pairs in the spinon condensation 
model Il5 20 , spin fluctuation in the nearly antiferro- 
magnetic ( AF) spin fermion model ^^, d-wave pairing 
fluctuations [^-^ , AF fluctuation- mediated pairing in- 
teractions in the d-p mode l |2q| , d-wave pairing fluctua- 
tions and AF fluctuations p6(|, and so forth. 

Among these candidates, the decrease of Tg (the tem- 
perature below which the pseudogap phenomena are ob- 
served) with hole doping seems to suggest a possibility 
that the pseudogap mainly originates from AF fluctu- 
ations at least at temperatures much higher than the 
superconducting transition temperature (Tc). The AF 
long-range-ordered phase occurs in the vicinity of the 
half-filling, where the real gap is open. It is likely to 
have a pseudogap structure in the DOS in its proximity 
due to strong AF fluctuations. In this mechanism, the 
decrease of Tg can be explained naturally, since the AF 
fluctuations decrease with doping. On the other hand, if 
we assume that the pseudogap is due to pairing fluctua- 
tions, temperature Tg can be regarded as the temperature 
at which pairing fluctuations begin to occur. Therefore, 
when Tg increases, it is natural to consider that Tc should 
increase as well. However, such behavior is inconsistent 
with the observed behavior of the opposite doping de- 
pendences of Tg and T^ in the underdoped region. 

In this paper we examine the pseudogap due to AF 
fluctuations. The reduction of the DOS near the Fermi 
surface should suppress Tc in the underdoped region, and 
thus Tc has a peak as a function of hole concentration. 
We intend to describe a minimum theory which repro- 
duces the phase diagram of a HTSC. Hence, we omit some 
of the details which do not change the qualitative behav- 
ior of Tc. For example, we adopt the static approximation 



for the spin fluctuations. The static approximation does 
not produce the imaginary part of the self-energy and the 
broadening of the one particle weight. However, the real 
part of the self-energy could reproduce the reduction of 
the DOS near the Fermi surface, which is the most in- 
tuitive deflnition of the pseudogap. Therefore, the static 
approximation is sufficient for our purpose. In addition, 
Schmalian et al. E^ argued that the characteristic fre- 
quency of the spin fluctuations LUgf is much smaller than 
the temperatures of interest for HTSCs. 

In our formulation, AF fluctuations are taken into ac- 
count through a renormalization effect. The importance 
of the renormalization effect in the doping dependence of 
the Tc of HTSC was discussed on the basis of the two- 
dimensional Hubbard model u% . It was shown that Tc is 
reduced considerably near the boundary of the AF long- 
range-ordered phase by the strong renormalization effect 
due to AF fluctuations. However, the pseudogap was not 
taken into account sufficiently. 

This work was extended to include the pseudogap due 
to AF fiuctuations in a quasi-one-dimensional (QID) 
Hubbard model as a model of QID organic superconduc- 
tors [^. It was shown that the pseudogap suppresses 
Tc markedly near the spin density wave boundary. As a 
result, phase diagrams of QID organic superconductors 
in the pressure and temperature plane were semiquanti- 
tatively reproduced. 

In HTSCs, however, the same approach based on one 
of the microscopic models is difficult. The interlayer cou- 
pling is much smaller and the temperature range of inter- 
est is much higher in HTSCs than it is in the organic su- 
perconductors. Hence, the thermal fluctuations are very 
strong in HTSCs. This situation makes a quantitative 
argument difficult. In addition, there is no common con- 
sensus as to which microscopic model is appropriate for 
HTSCs: a single-band Hubbard, d-p, t-J models, and 
so forth. Therefore, we treat phenomenologically de- 
termined AF fluctuations from the experiments instead 
of calculating them microscopically, and concentrate on 
their qualitative features. 

In the calculation of Tc, for simplicity, we assume that 



the coupling constant of the pairing interactions does 
not depend on the doping. Such interactions may be 
attributed to those mediated by phonons. It may ap- 
pear that the d-wave pairing does not occur with phonon- 
mediated interactions. However, it is shown that pairing 
interactions mediated by screened phonons can give rise 
to a d-wave pairing superconductivity in the presence of 
AF fluctuations p9| . Here, we do not specify the origin 
of the pairing interactions. To some extent there may be 
a contribution from the exchange of AF fluctuations 

We calculate the electron self-energy up to the one-loop 
approximation as 

E.(fc) = -T^ N-' ^ Yl ^-' (^' ^')G.'(fc'), (1) 

n' k' cr' 

with k = (k, iwn), where V^a' is the effective interaction 
due to the exchange of the fluctuations and Ga{k') is the 
renormalized electron Green's function. 



G,(fc) 



I 



Jn -Ck - Scr(fc) +^' 



(2) 



We consider a two-dimensional tight-binding model with 
the electron dispersion 



Ek — — 2t(cos kx + cos ky) — At' cos k^ cos k, 



1" 



(3) 



where t and t' are the nearest- and next-nearest-neighbor 
hoppings, respectively. 

We express the effective interactions due to exchange 
of magnetic fluctuations as 



nk,k') 



Vo 



(k - k' - q„ 



)^+%^' 



(4) 



within the static approximation. Here, Vq and qo are the 
phenomenological parameters. Since the effective inter- 
actions due to the exchange of magnetic fluctuations are 
proportional to the spin susceptibility x(k — k'), they 
must have a sharp peak at k = k' -I- q,„, where q^'s are 
the nesting vectors near (tt, tt) for the AF fluctuations 
which give the largest value of x(q). 

There are four nesting vectors such as q^ = (tt ± 5, tt) 
or qm = (tt, TT ± (5), when the incommensurate AF fluc- 
tuations occur iQ. In eq. (^, we take a nesting vector 
such as TT > qmx > Qmy > for kx — k'x > ky — k'y > 0. 
For other regions of k^ — k'^ and ky — k'y, q^ is chosen so 
that T^(k, k') satisfies the symmetry condition. The AF 
correlation length ^ = 1/qo diverges at the critical hole 
concentration n^. We take 
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7r(Mc-A^)i/2' 



(5) 



which diverges at a critical chemical potential Hc- For 
small hole concentrations, the chemical potential n is 
roughly proportional to the hole concentration nh- Thus, 
the above form roughly implies ^ ex 1/ ^/n^, which is plau- 
sible in the sense that the AF correlation length is di- 
rectly related to the average distance between holes [pi| . 



Thus, we finally obtain an equation to solve in a com- 
pact form 



1 



r-l 



S,(k) = -iV ^ 
k' 



F(k,k')tanh( 



Ck' -I- Ser(k') - ^ 
2T 



)■ (6) 



With the static approximation of eq. (^, we only need 
to solve the real part of the self-energy Thus, it is found 
that the pseudogap appears as shown below, for example, 
in the DOS, which is given by 



Pi')= /^^^(^-^~k), 



(7) 



where we have put ek = fk + ^^(k'). 

We solve the self-consistent equation for large discrete 
512 X 512 momentum points in the first Brillouin zone. 
We confirmed that the results of 512 x 512 practically 
coincide with those of 256 x 256, within the width of 
the Hues in our figures of DOS and T^. Thus, practi- 
cally, the system size 512 x 512 can be regarded as being 
within the thermodynamic limit at the present tempera- 
ture. We interpolate the obtained self-energy linearly in 
the momentum space so that the first Brillouin zone has 
4096 X 4096 points for the calculation of DOS. 

We consider the case of nearest-neighbor hoppings {t ^ 
and t' = 0). We choose Vq = -0.4 and T = 0.01, and 
take C — 30 and fic = —0.05 in eq. (||). We use the units 
in which t = 1. The value of /ic gives n^ « 0.02 — 2%, 
where n^ denotes the critical hole concentration of the 
AF long-range-order. 
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FIG. 1. Doping dependence of the AF correlation length 
^ (triangle), and the resulting densities of states psi^p) and 
Pdi^p), (circles and squares, respectively). The solid lines are 
guides for the eyes. Pdi^p) is the effective DOS for the d-wave 
pairing which is averaged with the weight [cosk^ — cosky)^. 

Figure |l] shows the n^ dependence of the AF correla- 
tion length ^, and the resulting DOS at the Fermi en- 
ergy. In this figure, pj^ is the effective DOS for the d- 
wave pairing in which the order parameter is A(k) ex 
(coskx — cos ky), 



Pd{^.) = 






fi) (cos kx — cos ky)'^ 



(8) 



(A(k) 



Ps is that for the s-wave pairing, 
is equivalent to the usual DOS. 

From the obtained DOS, we calculate T^ 
ventional weak coupling formula 



Tc = 1.13 uju e 



-l/Ao 



const.), which 
by the con- 

(9) 



with Aq = gpaiiJ,)- We choose a value of the cutoff en- 
ergy of bosons which mediate the pairing interactions as 
wd = 1000 K. The value of g is chosen so that Tc takes 
reasonable values such as ~ 40 K or '--' 90 K at the peak. 

Figures g and |^ show the doping dependence of the Tc 
for the s-wave pairing and d-wave pairing, respectively. 
From the sensitivity of the singular exponential form of 
eq. (P), we obtain marked suppressions of Tc near the AF 
boundary (rih ^ ^^h) ^^^ thus peak structures around 
TT-h ~ 0.1 = 10%. In particular, for the d-wave pairing, 
the peak is narrow, and seems to coincide with the ex- 
perimental phase diagram of HTSC. 

In order to be more precise, we should take into ac- 
count the temperature dependence of the DOS in the es- 
timation of Tc. However, for nh ^ n'^, the AF correlation 
length ^ does not strongly depend on temperature at low 
temperatures [0. When its temperature dependence is 
ignored, the DOS are almost independent of the temper- 
ature. For example, it can be confirmed that the results 
for T = 0.005 are almost the same as those for T = 0.01. 
On the other hand, even when the temperature depen- 
dence of ^ is taken into account, it does not change the 
qualitative result. It only emphasizes the peak of the Tc, 
because the pseudogap becomes deeper for longer ^ and 
it reduces Tc more. 
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FIG. 2. Doping dependence of Tc for the s-wave pair- 
ing. Circles, triangles and squares indicate the results for 
g = 2.3, 1.7 and 1, respectively. Solid lines are guides for the 
eyes. 




FIG. 3. Doping dependence of Tc for the d-wave pair- 
ing. Circles, triangles and squares indicate the results for 
g — 0.9, 0.7 and 0.5, respectively. Solid lines are guides for 
the eyes. 

Figure |j shows the DOS of the underdoped region 
(nh « 0.0237) and that of the overdoped region (nh « 
0.227). It is seen that the pseudogap is deep and clear 
in the underdoped region, while it becomes very shallow 
in the overdoped region. We also find that the pseudo- 
gap becomes large in the direction of (±7r, 0) and (0, ±7r), 
while it becomes small in the direction of (±7r, ±7r) 133]. 

It is straightforward to extend the present calculation 
to the t' ^ case. When we assume t' — —0.2 and a 
similar doping dependence of the AF correlation length 
which decreases with doping, we also have a peak struc- 
ture of the Tc near rih ~ 0.13 |32|. This peak struc- 
ture is also consistent with the experimental phase dia- 
grams except that the peak becomes very steep due to 
the van Hove singularity. However, it is easily verified 
by a phenomenological consideration that such steepness 
of the peak is not at all essential. If we assume that a 
large t' is appropriate for the HTSC, the Fermi surface 
crosses the van Hove singularities at a large hole con- 
centration. Thus, the DOS at the Fermi energy and Tc 
have a sharp peak around the hole concentration. How- 
ever, such singular behavior or marked enhancement of 
the DOS with hole doping has not been observed in any 
experiments, for example, the specific heat and suscepti- 
bility measurements. This suggests that the singularity is 
eliminated finally by some many-body effect. If we take 
this many-body effect into our phenomenological model 
beforehand, models with t' « might be reasonable in 
order to describe a realistic situation, because it is con- 
sistent with the experimental facts that no singularity 
occurs in the overdoped region, and the AF correlation 
length increases near the half-filling. 
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FIG. 4. DOS of the underdoped region /i = —0.06 and 
njj ~ 0.0237 (solid line), and that of the overdoped region 
/i = -0.6 and < « -0.227 (broken line). The dotted line 
indicates the DOS when Vo = 0. The vertical thin solid and 
broken lines indicate the chemical potentials /i = —0.06 and 
fj, — —0.6, respectively. 

We omit vertex corrections and dynamical effects. 
However, since these effects should reduce Tc in the un- 
derdoped region, it is likely that they only emphasize the 
peak structure of Tc. Also, the strong coupling effects 
would not be essential regarding the qualitative behavior 
of the Tc. The essential point is that Tc is very sensi- 
tive to the DOS, as explicitly written in the singular ex- 
ponential form eq. (|9|) within the weak coupling theory. 
The quantitative improvement of the theory by including 
these effects remains for future study. 

In conclusion, we propose a mechanism which explains 
the peak structure in the doping dependence of Tc of 
HTSCs. It is shown that the pseudogap due to AF fluctu- 
ations suppresses Tc in the underdoped region and even- 
tually destroys the superconductivity at a finite doping. 
On the other hand, the pseudogap phenomena is less pro- 
nounced near the optimum doping at which Tc is maxi- 
mum. Then, DOS is large since the Fermi energy is at 
the shoulder of the van Hove singularity. This leads to 
high Tc. The decrease of Tc in the overdoped region is 
due to the decrease of the DOS [pSipipgfl , since the Fermi 
energy is distant from the van Hove singularity. 

The authors would like to thank Professor D. Rainer 
and Professor J. Friedel for useful discussions and encour- 
agements. This work was partially supported by a grant 
for CREST from JST. 



[2] 

[s; 
[4; 

[5] 
[6] 

[7 

[9 

[10 

[11 
[12 

[13 

[14 



[is; 
[16; 

[17 



[is; 

[19 

[20 
[21 



[22' 
[23^ 

[24; 

[25; 

[26' 
[27 

[28] 



[29] 
[30] 
[31] 

[32] 
[33] 
[34] 



H. Ding et al: Nature (London) 382 (1996) 51; H. Ding 
et al: Phys. Rev. Lett. 78 (1997) 2628. 
J. M. Harris et al: Phys. Rev. B 54 (1996) R15665. 
T. Sato, et al: Phys. Rev. Lett. 83 (1999) 2254. 
Ch. Renner et al: Phys. Rev. Lett. 80 (1998) 149. 
N. Miyakawa et al: Phys. Rev. Lett. 80 (1998) 157; ihid 
83 (1999) 1018. 

T. Ekino et al: Phys. Rev. B 60 (1999) 6916. 
J. W. Loram et al: Phys. Rev. Lett. 71 (1993) 1740. 
H. Yasuoka et al: Strong Correlation and Superconduc- 
tivity (Springer- Verlag, Berlin, 1989), p. 254. 
W. W. Warren et al: Phys. Rev. Lett. 62 (1989) 1193. 
M. Takigawa et al: Phys. Rev. B43 (1991) 247. 
H. AUoul et al: Phys. Rev. Lett. 70 (1993) 1071. 
K. Ishida et al: Phys. Rev. B58 (1998) R5960. 
Neutron results are reviewed by P. Bourges, in The Gap 
Symmetry and Fluctuations in High Tc Superconductors, 
ed. J. Bok, et al (Plenum, New York, 1998), p. 349. 
Y. Suzumura, Y. Hasegawa and H. Fukuyama: J. Phys. 
Soc. Jpn. 57 (1988) 2768. 

G. Kotliar and J. Liu: Phys. Rev. B 38 (1988) 5142. 
T. Tanamoto, H. Kohno and H. Fukuyama: J. Phys. Soc. 
Jpn. 67 (1994) 2739; H. Fukuyama, H. Kohno and T. 
Tanamoto: J. Low Temp. Phys. 95 (1994) 309. 
P. Lee and X.-G. Wen: Phys. Rev. Lett. 78 (1997) 4111. 
M. Trivedi and M. Randeria: Phys. Rev. Lett 75 (1995) 
312. 

V. Emery and S. Kivelson: Nature 374 (1995) 434. 
J. Schmalian, D. Pines and B. Stojkovic: Phys. Rev. Lett. 
80 (1998) 3839; A.V. Chubukov and J. Schmalian: Phys. 
Rev. B 57 (1998) R11085. 

J.R. Engelbrecht et al: Phys. Rev. B 57 (1998) 13406. 
T. Ichinomiya and K. Yamada: J. Phys. Soc. Jpn. 68 
(1998) 981. 

T. Hotta, M. Mayr and E. Dagotto: to be published in 
Phys. Rev. B. 
A. Kobayashi, A. Tsuruta, T. Matsuura and Y. Kuroda: 

J. Phys. Soc. Jpn. 68 (19 99) 2506. 

S. Onoda and M. Imada: |cond-mat/9908195 . 
H. Shimahara and S. Takada: J. Phys. Soc. Jpn. 57 
(1988) 1044. 

H. Shimahara: J. Phys. Soc. Jpn. 58 (1989) 1735; 
H. Shimahara: Proceeding of the Physics and Chem- 
istry of Organic Superconductors, eds. G. Saito and S. 
Kagoshima (Springer- Verlag, Berlin, Heidelberg, New 
York, 1990), p. 73. 
J. Friedel and M. Kohmoto: |cond-mat/99Q1065t I. 



[1] D. S. MarshaU et al: Phys. Rev. Lett. 76 (1996) 4841; 
A. G. Loeser et al: Science 273 (1996) 325. 



Chang, J. Friedel and M. Kohmoto: |Eond-mat/9908214 . 
Y. Hasegawa and H. Fukuyama: Jpn. J. Appl. Phys. 26 
(1987) L322. 

R. J. Birgenneau and G. Shirane: Physical Properties of 
High Temperature Superconductors I, ed. D. M. Ginzberg 
(World Scientific, 1989). 

H. Shimahara, M. Kohmoto and Y. Hasegawa: unpub- 
lished. 

H. Fukuyama and Y. Hasegawa: J. Phys. Soc. Jpn. 56 
(1987) 1312. 

H. Shimahara and S. Takada: Jpn. J. Appl. Phys. 26 
(1987) L1674. 



